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ABSTRACT 


An exaalnatlon of the deternlniatic predictebllity for tropics and middle- 
latitudes separately Indicates that the theoretical upper limit of deterministic 
predictability for low latitudes is shorter than that for middle latitudes. 

Most of the day-to-day fluctuations in the tropics are determined by the growth 
and decay of condensation driven instabilities for which the amplitudes get 
equilibrated rapidly. The errors of observations are already closer to the 
mixiauM possible error for useful predictability. Therefore it takes only a 
few days for an initial error to grow to a magnitude comparable to the climato- 
logical variance. 

Variability of time averages in low latitudes is mainly determined by the 
location and intensity of the large-scale Hadley and Walker circulations. 

Since these are largely Influence by the slowly varying boundary conditions of 
sea surface temperature and soil moisture, and since synoptic instabilities 
are not strong enough to change drastically the large scale flow, therefore, 
there is larger potential for predictability of monthly and seasonal means in 
low latitudes. Since the tropical heat sources can also influence the varia- 
bility of the middle latitude circulations, under favorable conditions, time 
averages for middle latitudes can also be potentially predictable due to their 
interaction with the tropical heat sources. 

I It is conjectured that for short and medium range deterministic prediction, 
a prescribed diabetic heating field due to moist convection may be more useful 
than their explicit calculation from the evolving flow. Inadequacies of the 
current paramaterisation techniques rapidly degrade Che motion field, which in 
turn produces more unrealistic heating fields giving rise to still more 
unrealistic flow patterns. 
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Introduction ; 

The theoretical upper Unit for deteministic prediction is nainly deter- 
nlned by the growth rate and equilibration of doninant instabilities which 
exist for a given observed state of the atnosphere* An uncertainty in the ini** 
tlal conditions grows with the characteristic grotrth rate of the fastest growing 
instabilities. Nonlinear interactions aaong different scales of motion help 
spread this instability to all the scales present in the flow. For a sinple 
hydrodynaalcal system (barotroplc fluid without 6 effect and without mountain ) , 
Lorens (1969) showed that diffe.cat scales ur notion have different ranges of 
predictability and that the theoretical upper limit of predictability ranges 
from a few days to a few weeks. Several general circulation model studies have 
been carried out (Chamey et al. , 1966; Snagorinsky, 1969) to determine the 
theoretical upper limit of predictability. In these studies a general circula- 
tion model is first integrated with observed (or model simulated) initial con- 
ditions. Initial conditions are then perturbed by superimposing a random field 
of meteorological variables, and the model is integrated again keeping every- 
thing (boundary conditions, physical patameterisatlons, model, etc.) identical 
to the first run. Random perturbations in the initial conditions are assumed 
to simulate uncertainty in the definition of the Initial stated. Departures 
beween the two model Integrations are studied as a function of time to determine 
the growth rate of errors. On the basis of the rate of growth of the globally 
averaged root mean square error, it has been suggested that the doubling tisw 
fur the error is about three to five days. It should be noted that this approach 
only gives an estimate of the doubling time, and in order to determine the upper 

^ The inadequacy of the present observing systems gives rise to, in addition 
to random errors, large systoatic errors over the date-void areas. 


limit of predictability, one auet aseuae a value for the Initial error and 
another value for the final tolerable error* 
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In thla paper, we have shown that the error growth rates and their equili- 
bration depend upon the dynaalcal reglae, the aagnlttide of the error, the Ini- 
tial condition, and the aeteorologlcal variable. Since the nature of dynasdcal 
Instabilities which re aost laportant for day-to-day fluctuations In the 
tropics are very different froa those In the middle latitudes. It Is not appro- 
priate to examine the growth rate of combined error. In fact, even In the 
middle latitudes. It Is reasonable to assuae that the growth rates In the 
regions of highly active stora-tracks will be larger than those In the regions 
of quasl-peraanent anticyclones. These arguments are valid only for Initial 
growth rates and subsequently error growth rates will also be determined by 
Interaction aaong different scales and different dynaalcal regimes. For example, 
•for a suitable structure of the xonal flow, middle latitude effects can propa- 
gate to the lower latitudes and vice-versa, and for time periods longer than 
the time scale of these Interactions, one must consider the combined effects. 

The tropical circulation Is dominated by quasl-statlonary heat sources and 
associated Hadley and Walker type circulations. The space and tlae scales of 
these circulations are auch larger than the space and tine scales of tropical 
dlaturbances (easterly waves, depressions, cyclones, etc.). Therefore It Is 
less likely that the large scale tropical clrculatlona can be made unpredictable 
by their Interaction with the saall scale tropical disturbances. Although 
there Is evidence to suggest that the Initiation of the tropical disturbances 
Is aalnly due to dynaalcal Instabilities (barotroplc or combined barotroplc- 
barocllnlc). It Is well known that the further growth and the maintenance of 
tropical disturbances Is aalnly due to the latent heat of condensation, and 
CISK Is the aost dominant dynaalcal mechanism to explain their energetics. 
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Th« naximum aaplituda of theie diaturbancM la aqullibratad by availability of 
■olatura and efficiency of utilixing the available aoiature (Stevena at al« , 
1977; Shufcla, 1978). Ponution, growth and aoveaMnt of theae diaturbancea 
ia very iaportant for ahort-range forecaating in the tropica. Condenaation 
driven inatabilitiea (CISK) have large growth ratea but their aaplitudea eciil- 
ibrate rather quickly. Ihia unique nature of tropical diaturbancea can be a 
■ajor obatacle for ahort-range prediction in the tropica. Since the apetial 
acale of theae diaturbancea ia only about 2000-3000 ka and eabedded cloud 
cluatera are even analler, it ia required to have a aufflciently high reaolution 
to define their atructure. 

The fluctuationa of large-acale Hadley and Walker circular iona arc affected 
by the alowly varying boundary condltiona of aea-aurface teiq>erature (SST) and 
aoil aoiature. It la therefore conceivable that the tlae averagea (at leaat 
for the tiae acalea over which SST and aoil aoiature do not change rapidly) of 
the tropical circulation aay be aore predictable. Changea in SST can change 
the location of aacendlng branchea of the Hadley and Walker cella and atructure 
of diabetic heating fielda and therefore produce large changea in the tiae 
averaged precipitation and aaaociated circulation. For longer tiae acalea 
one ixuat alao conaider the Interaction between tropical and alddle latitude 
clrculationa. Fluctuationa in aoaentua and heat fluxea due to aiddle latitude 
eddy activity can affect the Intenaity of Hadley circulation and if the fomer 
are unpredictable » the latter would alao be unpredictable. Baaed on the above 
conalderationa , Chamey and Shukla (1981) auggeated that the tiae-averaged 
■onaoon clrculationa are aore predictable than the tine-averaged aiddle latitude 
circulation. 

The aiddle latitude circulation ia doalnated by atationary planetary 
wave a forced by orography and diabatic heat aourcea, transient long wnvea and 
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baroclinically unstabla synoptic seals wavss. Fast growing barocllnlc Insta- 
bllltlsa Intsract with ths largsr scalss and aaks than unprsdlctabls. Prsdlc- 
tlon of tlM avsragss (vis. for nnnthly asans) for alddls latltudss, thsrsfors. 
Is aors difficult than for low latitudes. Unllks ths tropical cass, ths Inter- 
action between barocllnlc eddies and planetary waves In the alddle-latltudes 
Is Buch stronger and therefore the potential for long range predictability Is 
saaller. 

However, the potential for short range deteralnlstlc prediction In alddle 
latitudes Is auch better than that for the low latitudes. The aaxlaua aaplltude 
of equilibration for alddle latitude eddies Is large enough so that the Initial 
uncertainty has to grow for several days before It becoaes coaparable to the 
aagnltude of the day-to-day fluctuations we wish to predict. 

In Part I of this paper we have utilised the 6LAS (Goddard Laboratory for 
Ataospherlc Sciences) cllaate aodel to exaalne the growth rate for Initial ran- 
doa errors for low and alddle latitudes separately. We have also coapared the 
errors due to random perturbation with the aaxlaua tolerable error (given by 
the standard deviation of the total tlae series) for each grid point and pre- 
sented the results of sonal averages. It Is found that for a reasonable value 
of the Initial error, the error becoaes about half of the aaxlaua error within 
2-3 days for the tropics and 5-7 days for the alddle latitudes. 

In Part II of this paper we have examined the effects of SST fluctuations 
between the equator and 30*N. He find support for the hypothesis proposed by 
Chamey and Shukla (1981) that the fluctuations of the bou'>4ary conditions 
can explain a significant part of the Innerannual variability of the tropical 
ataosphera. We have ezaalned the effects of SST fluctuations between the 
equator and 30*11. It Is found that the Int^ranaual variability of SST can 
axplaltt a significant part of tha variability of tha tropical ataosphera. 
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Id Part 111 of this paper we have sumariaed the reaulta of nv*;aeTical 
experlMnts carried out with global general circulation wodela which auggeat 
the Importance of boundary conditions for predictability of the time-averaged 
tropical circulation. 

Finally, we have conjectured that It la the quaal-ateady componenta of the 
tropical heat aourcea which are Important In affecting the circulation over 
middle latitudea, and therefore, although the tropical atmoaphere Itaelf la 
determlnlatlcally unpredictable for abort and medium range, Ita Influencea on 
middle latitudea could be calculated by preacrlblng the obaerved atructure and 
Intenalty of the tropical heat aourcea. 
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Rundoa Ttor, Cll— tologlcal grror mad PTti«t«nct Ttor, 

For a tiaa aarlas , tha varlanca (o^) la dafinad aa 

- 1 I (x^ - x)^ whara x - N l*t 

jm 1 ' 

If la Cha naan aquara dlffaranca batwaan all poaalbla paira of x^, it can 
ba ahown that 

E - /To 

E la a atatiatlcal aaaaura of arror betaaan two randomly choaan valuaa and a la 
a atatiatlcal maaaura of arror If tha long tarm aaan (climatology) waa aaatnad 
to bo tha foracaat for aach day. Thaaa daflnltlona ramaln valid both for a 
point (atatlon) and for a apatlal domain. The paraiatanca arror at any point 
changaa with time and tha magnitude of tha arror dapanda upon day-to-day atmo- 
apharlc fluctuatlona. E at any point, for aufflclently large aamplaa, can be 
a atatiatlcal maaaura of paraiatanca arror bacauaa It than bacomaa equivalent 
to ualng all poaalbla paira. Alternatively, If tha paraiatanca arror la cal- 
culated over a apatlal domain with a large number of polnta— large enough ao 
that tha apatlal dlatrlbutlon of arror at any time la comparable to tha time 
dlatrlbuf.lon of error of all tha polnta — tha maximum value of paraiatanca arror 
can be aatlmatad by E. Tha paraiatanca arror growa rapidly for a few daya 
before It attalna Ita maximum value. In tha abaanca of a long time aarlaa, tha 
maximum value of paraiatanca arror can ba uaad to aatlmata tha climatological 
atandard deviation. 

Tha theoretical upper limit of predictability la generally conaldarad to 
ba tlia time taken for tha arror to ba comparable to the error between two 
randomly choaan modal atataa. Since thla arror la /T tlmaa larger than tha 
atandard deviation of a climatological foracaat, which doaa not require any 
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•kill in prediction, it it too large an error to be a Maeure of ueeful predic- 
tion. In thie paper, «e ehall aeeuae o to be the aaxiBua tolerable error for 
calculating the lisita of predictability. It ia known that the atandard devla- 
tiona of day-to-day fluctuationa have a well defined geographical atructure, 
in particular, o at low latitudea ia analler than that at aiddle and high lati- 
tudea. It ia therefore aore appropriate to eatiaate predictability by coaparing 
the aagnitude of prediction error and o at each point aeparately. 

The upper liait of predictability ia deterained by ttM relative aagnitudea 
of the error growth rate and the aaxiaua day-to-day variability. If the growth 
rate ia aaall , it takea longer for the error to be cMiparable to the aaxi>' la 
error; if variancea due to day-to-day fluctuationa are aaall, it t^Ma only a 
few daya for the error to be coaparable to the aaxiaua error. However, the 
error growth rate itaelf ia not conatant with tiae. Saall errora grow aore 
rapidly than the large errora. Therefore, a aaaller growth rate doea not necea- 
•arily iaply longer predictability becaua^ aaall growth rato can be alaply due 
to the fact that the error la already cloae to ita aaxiaua value. 
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Part I» Short rang* pfdlcf btl^.ty of ttw tropical ataotphT*. 


Wa hava axamlnad tha arror growth aaK>ng four ganaral circulation oK>dal 
runa for which tha initial conditiona wara randomly parturbad. Tlia aodr.l waa 
firat intagratad for 45 daya atarting froa tha Initial conditiona of aiddla of 
Juna. Thraa additlonil runa wara aada by parturbing tha initial conditiona of 
aaa laval praaaura, U» W, and T at aach of tha nina lavala of tha modal. Bach 
grid point waa randomly parturbad corraaponding to a gauaaian diatribution 
with aaro maan and atandard daviation of 3 m/a for u and v, 1*C for T and 1 mb 
for aaa laval praaaura. It ahould ba noted that thaaa arrora of praaant obaarv- 
Ing ayatam are comparable to day-to-day fluctuationa of tha tropical flow. Tha 
boundary conditiona ware identical for all tha four modal intagrationa . Vfa 
have calculated tha atandard daviation among the four runa at aach grid point 
for aach day. Theae modal runa ware carried out by Charnay at al. (1977) and 
daalgnatad aa predictability runa in thair paper. 

If any variable R for modal run a at grid point l,j at time t ia denoted 
•a 8^,j,t,m ** calculate tha error (a) among the four runs as follows: 

“ _ , 1/2 
•l.J.t - («ll.J.t.m - 1l,J,t)V(»‘-l) . (M-4) 


_ M 

where (R) - 1/M ^ (R) is tha average for tha four runa. 

wl have also calculated the standard deviation (a) at aach grid point among all 
tha modal states raallsad during 45-day integration of tha four runs. 

• i* a maasura of tha error of a climatological foracast. Ua would consider 
tha upper limit of predictability to ba tha time taken for a to ba comparable 


to a 
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Figures U and lb show the evolution of the errors for ses level pressure 
and U St 175 mb, respectively. In esch figure, the solid line, dsshcd line and 
dotted line refer to the error average over the 10* latitude belts centered at 
6*N, 30*N and 58*N, respectively. It Is seen that the error in the sonal velo- 
city at the end of one day is largest at 6*M, followed by 30*M and 58*N. The 
eras was true for the aerldlcnal velocity at 835 mb (not shown). The final 
values of error for U at 175 mb at 6*N, 30*H and 58*N are comparable because 
during the susHMr the tonal velocity and Its fluctuations at this level are 
very large for all the three latitudes; at 6*N due to an easterly Jet stream 
and at 30*N due to a subtropical jet stream. The final values of error for V 
at 835 mb (not shown) were comparable for 6*N and 58*N. This Is because the 
middle latitude summer circulation Is not as vigorous as the winter circulation 
and the day-today fluctuations In meridional velocity are not ao strong. The 
final value of error for sea-level pressure Is smallest at 6*N and largest at 
58*N. This Is due to somII day-to-day fluctuations of sea level pressure In 
the tropics and larger day-to-day fluctuations In the middle latitudes. 

Figures 2a, 2b, and 2c show the zonally averaged values tne ratio (e/s) 
for the Bonel velocity (U) at 175 nb, zonal velocity (U) at 835 mb and sea 
level pressure respectively. In each figure the zonal average of the ratio 
(e/s) Is shown for day 1, day 2, day 3 and day 7. The following conclusions 
for summer emerge from these figures. 

1) The error In the tropics becomes half of the climatological standard 
deviation within 2-3 days, wtiereas, for middle-latitudes. It takes 
abouc 5-7 days. 

11) The error growth rate for sea level pressure Is different from that 

for wind. In the tropics, the error In sea level pressure reaches Its 
maximum value much faster than that In the winds. The seme Is true 
for the temperature (not shown). This, of course, also depends upon 
the magnitude of the Initial error and the magnitude of s %rhlch depends 
upon the time variability of the parameter under consideration. 
Fluctuations of pressure and temperature are rather small In the tropics. 
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Saaller values of s in the tropics suggests that the day-to-day fluctua- 
tions are noc as large as in the niddle-latitudes. In the tropics, the errors 
of observation are already closer to the naximum permissible error for useful 
predictability and therefore it takes only a few more days for the error to 
be equal to the maximum error. 

Figure 3 shows the zonally averaged values of doubling time for different 
latitudes. The doublng time for the low latitudes is about 4-10 days which is 
higher than that for the middle and high laitudes. In particular, sea level 
pressure shows the largest doubling time in the tropics. This does not mean 
that the tropics have longer predictability because, as discussed earlier, the 
doubling tioKS are larger because the errors are already close to their maximum 
value. It can be seen from Figures 2a, 2b, 2c that the predictability is 
smallest for sea level pressure compared to U at 175 mb and 835 mb. 

Ue have repeated these calculations for the version of the GLAS model 
described by Halem et al. (1980). Starting from the initial conditions in the 
middle of June, the GLAS model was first integrated for 90 days. A second 
integration for 90 days was also carried out in which only Initial conditions 
of U and V fields at all the nine levels of the model were randomly perturbed. 
Each grid point was randomly perturbed corresponding to a Gaussian distribution 
with zero mean and standard deviation of 3 m/sec for U and V fields. Figures 
4a, 4b, 4c and 4d show the evolution of the initial error lot sea level pres- 
sure, U, V and T at 500 mb, respectively. In each figure, the solid line, 
dashed line and dotted line refer to the error averaged over a 10* latitude 
belt centered at 6*N, 30*N and 58*N, respecr.ively. It should be noted that 
although there was no initial error in the temperature and pressure field, the 
error in these fields for the first 4-5 days is the largest at 6*N. The same 
is true for errors in U and V components. This suggests that the rate of 


ll 


growth of Initial error is large in the tropical latitudes compared to the 
middle latitudes. The final asymptotic value of the errors in sea level pres- 
sure and 500 mb temperature field is the largest for 58°N and the smallest for 
6°N. This is due to large day-to-day fluctuations in sea level pressure and 
temperature in the middle laitudes compared to the low latitudes. Due to the 
smallness of coriolis parameter, large changes in pressure and temperature 
cannot be sustained in the tropics. The large values of the final error in U 
and V components are due to the subtropical jet stream near 30**N and tropical 
jet stream near during the northern summer. 

We have also calculated the error growth rate and predictability of the 
GLAS model for winter initial conditions. Figure 5a shows the zonally averaged 
root mean square error for geopotential height for nine pairs of model integra- 
tions in which the winter Initial conditions of U and V at nine levels of the 
GLAS model were randomly perturbed. It is seen that the error in the tropics 
is the smallest and the error in the northern hemispheric middle and high lati- 
tudes is the largest. The error in the southern hemisphere (summer) middle 
latitudes is larger than the tropics but not as large as the northern hemisphere 
(winter) middle latitudes. Figure 5b shows the zonal average of the ratio of 
the error and standard deviation of day-to-day fluctuations. It is seen that 
the ratio is large for the low latitudes compared to the middle latitudes. 

This indicates that although the error is smaller for the low latitudes, the 
magnitude of the day-to-day fluctuations is so small that it takes only a 
few days for the error to be about half of the standard deviation. 

The above results suggest that the short term predictability for the 
tropical atmosphere is limited only to a few days because of two reasons: the 
rate of growth of the initial error is large and the maximum possible value of 
error is small. If the initial errors (due to observational errors, lack of 
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observations or deficiency of analysis and Intlallsatlon schemes) were compar- 
able for Che low and Che middle laClCudes, Che limit of predictability for the 
low latitudes will be lower even If the growth rates were comparable. There 
is reason to believe that, since the low latitude Instabilities are driven 
mainly by condensation, they would be relatively more unpredictable because of 
the limitations on the physical parameterlsatlon of moist convective processes. 
This, we believe. Is especially true for the northern summer when moist convec- 
tion Is very Important In determining Che nature of day-to-day fluctuations. 

This would suggest that the results of such predictability studies, and there- 
fore also the results of actual numerical prediction, will depend upon parameter- 
lsatlon of physical processes of moist convection, boundary layer processes and 
cloud-radlatlon Interactions. The spatial scale of the tropical easterly waves 
and depressions Is also smaller than chat of mid-latitude barocllnlc eddies and 
thus they reqlure higher resolution for defining their structure. Even higher 
resolution would be required to define the mesoscale structures embedded In the 
tropical disturbances. On Che other hand, for middle latitude disturbances, a 
quantitative treatment of the main energetic processes of conversion from avail- 
able potential energy to kinetic energy Is more reliable because It depends 
upon the larv^'-'cale vertical velocity and temperature fields. 

Based on these considerations alone. It can be Inferred that even for an 
Idealized case of a uniform and high density of observations over the globe, 
the upper limit of deterministic predictability %rlll be shorter for the tropics 
than for the middle latitudes. The reality of the situation Is, of course, 
much worse. Tropical areas, even In the northern hemisphere, have far less 
upper air soundings compared to the northern hemispheric middle latitudes. 
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An additional Important problaa for abort ranga pradlctablllty la ralatad 
to tha traatnant of aountalna. Inadaquata traatnant of larga-acala orographic 
affacts (atatlonary forcing on a planatary acala) and aoall acala affacta (lea 
cyclogeneala, flow over and around mountains) Is one of tha serious daflclaucles 
of present NWF models. Although these effects are Important both at low lati- 
tudes and middle latitudes, It Is more difficult to treat the Interaction of 
local topography and tropical disturbances which are primarily driven by conden- 
sation. 
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P«t II. Predictability of monthly and seasonal means for tropics. 

The interannual variability of monthly and seasonal iMans is determined 
by the combined effects of the internal atmospheric dynamics and the slowly 
varying 'external' boundary conditions of sea surface temperature (SST), soil 
moisture, snow/sea ice, etc. Although solar heating is the only truly external 
forcing to the atmosphere, the boundary conditions at the earth-atmosphere 
interface can be assumed to be external at least for those time scales for which 
they change slowly compared to the atmospheric dynamics. There is observational 
evidence to suggest that these boundary conditions in fact do change slowly 
enough so that they may be thought of as constant in time up to even a season, 
which is a longer time scale compared to the time scale of synoptic scale insta- 
bilities which are most important for day-to-day fluctuations. In the absence 
of any changes in these boundary conditions, the internal dynamical changes 
(which include the earth's topography and the land-sea contrast) will produce 
interannual variability of monthly and seasonal means. If the total interannual 
variability could be explained by the internal dynamics alone, the potential 
for predictability of time averages would be rather low. It is difficult to 
isolate the contribution of these two factors because the boundary conditions 
themselves are affected by the atmospheric dynamics, and observed variability 
is due to complex interaction between the two (Straus and Halem, 1981). Real- 
istic physical models of the atuosperlc circulation can be useful in conducting 
controlled numerical experiments to determine the relative roles of internal 
dynamics and boundary conditions. 

Chamey and Shukla (1981) have suggested that the time-averaged monsoon 
circulation is potentially more predictable than the middle latitude circulation. 
This is so because the large-scale monsoon circulation lb stable with respect 
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to dyn«alc instabilities which develop in the aonsoon flow, and fluctuations 
in the boundary conditions hav" significant effects on the tlae-averaged mon- 
soon flow. This conclusion was arrived at by exaalning the variability aaong 
the monthly Man (July) circulations of four model runs for which the boundary 
conditions were kept identical but the Initial conditions were randomly per- 
turbed. It was found that, although the observed and the model variabilities 
were comparable for middle and high latitudes, the variability among the four 
model runs for the monsoon regions was far less than the observed interannual 
variability of the atmosphere. This led us to conclude that the remaining 
variability could be due to the boundary conditions. This was also consistent 
with the results of several observational studies and numerical experiments 
which showed that the changes in SST or soil Misture at low latitude produce 
significant changes in the atmospheric circulation. 

In this paper we have extended the work of Charney and Shukla (1981) by 
comparing the model variability due to internal dynamics and due to changes in 
the boundary condition of SST. One of the limitations of the earlier study 
vas the comparison of the model variability with the observed variability. 

While this must be the ultimate goal, it is more appropriate first to Intercom- 
pare two different properties of the saM model so that any deficiencies of 
the model Itself do not bias the conclusions. 

We have carried out a 45-day integration of the GIAS model starting from 
the observed initial conditions in the middle of June, and climatological mean 
boundary conditions of SST. We refer to this integration as control run (c). 

For the identical boundary conditions we have carried out three additional 
integrations for 45 days each by randomly changing the initial conditions of u 
and V at each of the nine levels of the model. The spatial structure of the 
random errors corresponded to a Gaussian distribution with sero mean and standard 
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deviation of 3 m/a for u and v separately. We refer to these three Integrations 
as predictability runs (Pi, ? 2 » ’Ps)* Although the statistical properties 

of the random errors mere the same for each predictability run, the actual grid 
point values were randomly different. We have also carried out three additional 
integrations for which, in addition to the randomly perturbed Initial conditions, 
the boundary conditions of SST between the equator and 30*N were replaced by 
the observed^ SST during July of 1972, 1973, 1974. We refer to these three 
integrations as boundary forcing runs (Bi, B£, and B 3 ). The differences between 
the climatological SST as used in the control and predictability runs, and the 
observed SST used in boundary forced runs is shown in Figures 6 a, 6 b, and 6 c. 
Although there are large systematic differences over a few grid points, the 
SST anomaly over most of the tropical oceans appears to be realistic. 

The variance (op)^ among C, Pp P 2 , and P 3 will give a measure of the 
natural variability of the model; the variance (og) «song C, Bp B 2 , and B^ 

will give a measure of the variability due to changes in the boundary conditions 

of tropical SST. We have also calculated the observed variances (o^)^ for ten 
years of observed monthly means. 

|7c - P)^ + (Pi ” P)2 + (P 2 - P)^ + (P 3 - P)^ i,j 

(aB)^l.j - 1/3 Qc - B)2 + (Bj - B)2 + (B 2 - B)^ + (B 3 - B)^^ l,j 

where P - (C + Pi + P£ + P 3>/4 

and 1 - (C + Bi + B 2 + B 3>/4 

where Cij , Pij , Bij denote the July mean at grid point ij . 

Figure 7 shows the plots of sonally averaged values of standard deviations 

Op> 03 > Oo ^ 1 ^ ratios Oo/op and Oq/ob* agreement with the results 

of Chamey and Shukla (1981), it is seen that the ratio oo/op more than 

3 The observed SST was kindly provided to us by Dr. E. Kraus. 
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two in the tropical latitudes and close to one in the Biddle latitudes. The 
new result is that the curve ob lin* nearly halfway between the curves oo 
Op. This suggests that about half of the remaining variability for this modal 
is accounted for by changes in sea-surface temperature between equator and 30*N. 

This supports our earlier hypottssis that the slowly varying boundary con- 
ditions play an important role in determining the interannual variability of 
time averages for the tropics. Additional effects of soil moisture or the Eurasian 
snow cover could possibly bring the oo OB curves still closer. It is 

however to be noted that the long period internal dynamical changes (troplcal- 
extratroplcal Interactions, etc.) also contribute to the interannual variability 
of time averages and it could never be possible to explain the total qq by 
boundary conditions alone. In one of the model integrations for 90 days des- 
cribed in the earlier section, it was found that large differences in the monthly 
means were produced by internal dynamics alone. 

We have also compared the model variability for the predictability and 
boundary forcing runs. Since the sea surface temperature anomalies for B^, B 2 , 
and B 3 have many common features, we have considered it more appropriate to 
calculate the changes in the monthly means due to boundary conditions (Eg) and 
due to random perturbations (Ep) as follows: 

Ej^ - 1/3 (C - (at each grid point i,J) 

k«l 

2 ^ 2 

and Ep - 1/3 ^ (C - (at each grid point 1,J) 

k -1 

Figure 8 shows the aonally averaged values of Eg and Ep for July mean geo- 
potential height at 300 mb. Curves for Ep and Eg are labeled as 'PREDICTABILITY* 
and 'SST AHOMALT* respectively. As is well known, the values of Eb nnd Ep are 
small for the low latitudes and large for the middle latitudes, however, the 
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ratio Eb/Ep is iK>re than two for low latitudes. Tha largast values of the 
ratio Eg/Bp occur between 20*N-20*S. Thare are secondary eaxioa near the poles 
which occupy very aaall surface are« snd they can be ignored. 

We have noted, but we have not explained, that although the SST anoaaly 
was iapoaed only betiieen the latitudes 0-30”N, the effects on circulation are 
seen in the southern hemispheric tropics also. This could be due to meridionally 
propagating Rossby waves forced by heating due to SST anomalies, and interhemi- 
spheric interactions associated with the fluctuations of Hadley cells. The 
results for geopotential height at 500 mb (cot shown) are very similar to the 
one shown in Figure 6 except that the peak at the equator is not as high. These 
results suggest that SST anomalies in the tropics produce large changes in the 
middle and upper troposphere. Changes in the organisation and intenr^ty of 
deep moist cnvection Introduce significant changes in the diabatic heating of 
the upper troposphere. For suitable structure and intensity of the prevailing 
notion fields, these effects can be further transmitted to the middle latitudes. 
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S«n»itlvity of th> tropical circulation to cluagaa In th« boundary conditions* 

We have carried out several experiaents to test the sensitivity of a glo- 
bal general circulation aodel to changes In the plowly varying boundary condi- 
tions of sea-surface teaperature, soil aoisture and albedo. It is found that 
SST anoaalies in the low latitudes produce significant changes in the low and 
the aiddle latitudes. For exaaple, it is found that the warn (cold) SST ano- 
aalies in the Arabian Sea increase (decrease) acnsoon rainfall over India and 
adjoining areas (Shukla, 1975). It is also found that a siaultaneous occurrence 
of warm SST anotsalies over the north equatorial Atlantic and cold SST anoaalies 
over the south equatorial Atlantic produce severe drought conditions over north- 
east Brasil (Hours and Shukla, 1981). A wans anoaaly in the north and a cold 
anomaly in the south generates a thermally driven circulation whose descending 
branch is over northeast Brasil and adjoining oceans. An SST anomaly over 
equatorial Atlantic also produces significant aiddle latitude response in the 
northern hemisphere. There are nus^rous other observati )nal and aodellng 
studies which have shown that the tropical SST anomalies influence the intensity 
of Hadley and Walker cells. Statistically significant relationships have been 
found between the southern oscillation, SST anosMlles in the equatorial Pacific 
and the upper air circulation over Northern Hemispheric islddle latitudes. This 
suggests that the tropical thermal forcings can contribute to the predictability 
of the aiddle latitude time averages, which otherwise are less predictable by 
themselves. 

We have also carried out sensitivity studies for ttro extreme conditions: 

In one case, the soil is dry (no evapotranapiration) , and in the other case, 
the soil is wet (potential evapotranapiration) over global land surfaces (Shukla 
and Mints, 1981). It is found that during northern summer, absence of land- 
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lurfacc cvapc transpiration ovar nnat of tha healaphara raducas tha aonthly 
oaan rainfall by about half. Tha only axcaptions ara tha aonsoonal raglons 
for which raduction in avaporation ovar land la aora than coapanaatad by tha 
incraasad aolstura f7'oc convarganca asaociatad with intansa hast Iowa which 
fora duo to intanaa hasting of tha ground and tha ovarlying air. In tha abaanca 
of land avaporation ao^it of tha aolar radiation is utilizad to hast tha ground 
diractly and tha ovarlying air is haatad by sansibla hasting. 

Thara ara also a faw obaarvational and nuaarical atudlas which hava shown 
that the extant and tha depth of the Eurasian snow cover during winter ia 
related to tha Intensity of tha Aaiatc aonaoon circulation during tha following 
suoaer (Hahn and Shukla, 1977; Yah at al. 1981). The actual aachaniaa is not 
quite clear but large snow can give riee to large soil moisture which will 
impede the heating of the ground needed for the oneat of the monsoons. 

These results collectively suggest that the boundary conditions at tha 
earth's surface may be a useful predictor for the tropical circulation and 
under favorable conditions they can contribute to the predictability of the 


middle latitudes 
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Part III. Influ«nc« of tropical forcing on thm circulation of th« •xtr«-tropic«l 
ataoaphara . 

(1) It la a fact of obaarvatlona that tha subtropical highs in tha North 
Atlantic and tha South Atlantic hava tha highast monthly aran saa laval prassura 
sinultanaously during the month of July. This is rathar ramarkabla bacausa it 
is contrary to «hat one wotild axpact from axpactad phasa lag for tha saasonal 
cycle in tha two hamispheras. A possible reason for simultaneous iutansifics" 
tion of North and South Atlantic highs may be that they arc forced by tropical 
heating which is also maximum during July. This argument can be ganaralisad to 
include all the subtrop/cal highs between 30*N and 30*S. Figures 9a and 9b 
show global naps of 16-;ear (1961-1976) naan saa laval prassura for January and 
July. Wa have astinatad tha intensity of five subtropical highs (North and 
South Atlantic, North and South Pacific, and South Indian Ocean subtropical 
highs) by naasuring tha area enclosed by 1020 nb isobar. Figure 10 shows the 
plots of tha intensity of the subtropical highs as, neasurad by tha total area 
covered by tha five subtropical highs and rainfa.i' *'atwean tlw latitude belt 
30*N and 30*S for January through Decanbar. Tha rainfall data is taken fron 
Jaeger (1976). It is seen that tha intensity of tha subtropical highs is 
closely related to tha anount of precipitation and therefore tha latent heat of 
condensation in the tropical belt. It is likely that the intarannual variabil- 
ity of the subtropical highs may also be related to tha Interannual variability 
of tha tropical precipitation. However, adequate precipitation data is not 
available to test this hypothesis. It should be pointed out that stronger sub- 
tropical highs would imply stronger trade winds and possibly stronger inter- 
tropical convergence; however, its relationship to rainfall would also depend 
upon the location of the subtropical highs. 
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(2) Figure 11 show* the latlcudinal distribution of the itatlonery vari- 
ance of geopotential height at 500 ab aa alaMilated by two vereiona of the GIAS 
cliaate nodal. The noat inportant difference between the two aiwilationa waa 
the change in evaporation and precipitation in the tropical belt. The evapora- 
tion and precipitation in the old GIAS clinate nodei. waa leaa than tie obaerved; 
the evaporation and precipitation in the w>dified GIAS clinate nodal la nore 
than *:l)Q old nodal and conparable to the obaervationa. The change in the ata- 
tionary variance in tlM niddle latitudea ia rather atriklng. Although precipi- 
tation and evaporation in the nid-latitudea in the nodified nodal were not 
Identicil to the old nodal, and therefore all the change cannot be attributed 
to the tropical heating, it la reaaonable to auggeat that the tropical heating 
ia cne of the inportant factora to influence the nodal ainulated atationary 
variance. 

(3) There ia aone obaervational evidence that peralatent deep noiat acti- 
vity over the tropica can influence the niddle latitude circulation within a 
few daya. Paegle (1981) hea ahown an aaaociatlon between atrong outflowa at 
200 ab between equator and 20* S, and increaaed tranalent kinetic energy at 
30*N during January and February, 1979. Bxanlnatlon of dally cloud plcturea 

and upper level aynoptlc weather nape leavea one with the inpreaaion (J. Wlnaton, 
peraonal coanunication) of atrong aaaociatlon between Intenae t.opical convec- 
tive activity and intenaification of aubtropical jet atrean. 

If tropical activity haa atrong influence on nld-latitude weather fluctua- 
tiona, and if tropical flow ia unpredictable beyond a few daya, doea thia ieply 
that ^t would alao lielt the ■J.d-latltude predictebility? It la our conjecture 
that lack of detereiniatic predictability for aynoptlc acalea in low latitudea 
doea not neceaarily ieply that tha large-acale tropical-ewtratropical intarac- 
tiona cannet be adequately aodelled. It ia the quaai-atationary tropical heat 
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tourcea which are tha aaln drivara of cha txopical-axtratroplcal talacoanactiona 
and produca larga axtratropical raaponaa. It la concalvaabla that by aultabla 
procaduraa of aaalnllatlng vartlcal valocity and diabatlc haatlng flalda darivad 
froB tha obaarvationa of cloudlnaaa and pracipltation, tha larga-acala tropical 
haat aourcaa can ba raaaonably dafinad. It ia alao conca^vabla that thair 
affacta can be accounted for by ptaacribing than for a few daya. 

(4) Kalnay-Rivaa and Halaa (peraonal cooiBunication) have carried out aav- 
aral foracaat axpariaanta during FGGE SOPI with the GLA5 foracaat aodal. For 
four of tha caaaa they inaartad tha obaarvad FGGE data every aix houra during 
tha foracaat between 20*S and 20*N. Thia waa intended to aiaulata a perfect 
foracaat in the tropica. They compared th* akill of the forecaata over North 
America for th,< normal forecaata and the forecaata with tropical FGGE data 
inaartad during the forecaata. In three out of four forecaata, inaartion of 
tropical data had little effect on the forecast akill over North America; how- 
ever, in one of the four caaaa there waa aubatantial improvement in tha foracaat 
for North America after tha first three daya of integration. Cloud picturea 
ahowed the preaence of an intenae tropical flux originating in the equatorial 
Weatern Pacific and penetrating deep into tha aid-latitudea reaching the mat 
coaat of America. 

Although more obaervational and numerical atudiea are needed to eatabliah 
the nature of tropical-mld-latitude interaction at auch 'faat' time acalea, 
preliminary reaulta auggeat that tropical phenomena, in aome aituations, may be 
important even for abort range NWP in middle latitudea. 

(5) There ia yat no quantitative aaaeaament of the relative roles of data 
aaaiallation and initialisation procedurea, parameterisation of moist convective 
heating, and poor definition of initial state for rapid degradation of tk*opi.:al 
forecasts. If moist convection parameterisation ia not realistic, it quickly 
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produces unrealistic flow which in tura leads to still nore unrealistic moist 
convective heating. This fast feedback between the motion field and the moist 
convective heating may be one of the most Important limiting factors for tropl' 
cal predictability. Large sensitivity of tropical forecasts to changes in the 
initial moisture field support this point. It is therefore conceivable that 
prescribing the condensation heating may be less damaging than explicitly cal- 
culating it. This could be true even for middle and high latitudes. 
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Conclusions ; 

Classical predictability studies using a general circulation model show 
that for short- and medium-range the tropical atmosphere Is less predictable 
than the middle latitudes. This Is primarily because of the following reasons: 

1) Tropical areas have less data and therefore more uncertainty In the Initial 
conditions. The error of observations Is comparable to the changes to be pre- 
dicted. 2) Host of the day-to-day fluctuations are due to condensation driven 
Instabilities which grow rapidly. The Initialization and assimilation techniques 
do not take Into account the role of diabetic forcing for tropical disturbances 
and It Is difficult to parameterize the physical processes of moist convection, 
boundary layer and cloud-radlatlon Interaction. 3) The standard deviation of 
day-to-day fluctuations Is relatively small because the amplitudes of the tropi- 
cal disturbances equilibrate rapidly. It takes only a few days for the Initial 
error to grow to be comparable to the standard deviation of dally fluctuations. 

On the other hand, the time averages (monthly and seasonal means) for the 
tropics have more potential predictability. This Is because they are largely 
determined by fluctuations In the slowly varying boundary conditions of sea 
surface temperature and soil moisture. Under favorable conditions they can 
contribute to the predictability of middle latitudes also. 

Evidence Is beginning to emerge that tropical heat sources can also Influ- 
ence the middle latitude circulation within a few days. If the tropical atmo- 
sphere were Intrinsically unpredictable for the short and medium ranges. Its 
Influence on the middle latitudes will also be unpredictable. We conjecture 
that It Is the quasi-steady component of the tropical heat source which affects 
the middle latitude circulation and therefore It should be possible. In principle, 
to prescribe It for a few days from the observations. This may require special 
techniques of Initialization and assimilation of tropical data of cloudiness 
(vertical velocity) and rainfall (heating). 
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Figures 


Fig. le. 

Fig. lb. 
Fig. 2a. 

Fig. 2b. 
Fig. 2c. 
Fig. 3. 

Fig. 4a. 

Fig. 4b. 
Fig. 4c. 
Fig. 4d. 
Fig. 5a. 

Fig. 5b. 

Fig. 6a. 

Fig. 6b. 
Fig. 6c. 


Root nean square error between four suaner model runs as a function 
of time for sonal velocity (m/s) at 175 mb. Solid line, dashed line 
and dotted line refer to an average over 10* latitude belt centered 
at 6*N, 30*H and 58”N respectively. 

Same as Fig. la but for sea level pressure (mb). 

Zonally averaged values of the ratio of root mean square error between 
four summer model runs and standard deviation of dally values for the 
same four runs for zonal velocity at 175 mb. Curve labeled as DAY 
1,2,3 ,7 refer to the ratio at the end of one, two, three and seven 
days. 

Same as Fig. 2a but for zonal velocity at 835 mb. 

Same as Fig. 2a but for sea level pressure. 

Doubling time (days) for the Initial error among the four predictability 
runs of sea level pressure (solid line), zonal velocity at 835 mb 
(dashed line) and zonal velocity at 175 mb (dotted line). 

Root mean square as a function of tiiM between a summer control and a 
predictability run for sea level pressure (mb). Solid line, dashed 
line and dotted line refer to an average over 10* latitude belt 
centered at 6*N, 30*N and 58*N respectively. 

Same as Fig. 4a but for zonal velocity (m/s) at 500 mb. 

Same as Fig. 4a but for meridional velocity (m/s) at 500 mb. 

Same as Fig. 4a but for temperature (*C) at 500 mb. 

Root mean square error between four winter model runs as a function 
of latitude and time (day) for geopotential height at 500 mb. 

Zonally averaged values of the ratio of root mean square bet%reen nine 
pairs of winter model runs and standard deviation of daily values for 
geopotential height at 500 mb. 

Difference between the climatological sea surface temperature used 
for the control run and observed sea surface temperature during July 
1972 between eq-^ator and 30*N. 

Same as Fig. 6a but for July 1973. 

Same as Fig. 6a but for July 1974. 
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Fig. 


Fig, 

Fig. 

Fig. 

Fig. 

Fig. 


7. Zonally averaged atandard deviation among monthly mean (July) aea 
level preaaure (mb) for 10 years of obaervatlona (Oq, thin solid line) 
dashed line) and four model runs «rlth Identical boundary conditions 
(Op, thin dotted line). Ihlck solid line atul thick dashed line 

show the ratio Oo/op end Oq/ob tespectlvely. 

8. Zonally averaged standard deviation for predictability runs (thin solid 
line), boundary forced runs (dashed line) and the ratio of boundary 
forced and predictability runs (thick solid line) for geopotentlal 
height at 300 mb. 

9a. Sixteen year (1961-76) mean sea level pressure (-1000 mb) for January. 

9b. Sixteen year (1961-76) mean sea level pressure (-1000 mb) for July. 

10. Variation of monthly mean rainfall (mm) averaged between 30”N and 
30*S (from Jaeger, 1976), and intensity of subtropical highs measured 
by number of (4* 1st. x 5* long.) grid points for which sea level 
pressure Is greater than 1020 mb (from Godbole and 'Jhukla, 1981). 

11. Stationary variance of January mean geopotentlal height at 500 mb 
simulated by two versions of the GLAS climate model, and observations. 




Fig. U. Root aean square error between four suaaer aodel runs as a function 
of tlae for zonal velocity (a/s) at 175 ab. Solid line, dashed line 
and dotted line refer to an average over 10* latlttide belt centered 
at 6*N, 30*N and S8*N respectively. 
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Fig. 2b. Saae as Fig. 2a but for zonal velocity at 835 ab 
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Fig. 6a. Difference between the climatological sea surface temperature used 

for the control run and observed sea surface temperature during July 
1972 between equator and 30*N. 
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fig* 6b. Same as Fig. 6a but for July 1973 
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Zonally averaged standard deviation for predictability runs (thin solid 
line), boundary forced runs (dashed line) and the ratio of boundary 
forced and predictability runs (thick solid line) for geopotential 
height at 300 nb. 
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Fig. 9«. Sixteen year (1961-76) mean 5ea level pressure (-1000 ab) for January 






Fig. 9b. Sixteen year (1961-76) aean sea level pressure (-1000 ab) for July 



zoo 


o 


o 

0 - 



{mj) n\fdNi¥d 
o 
o 



(ewoaoi sih/tod aidS> jo ddst^nn fc/7sr) 


MONTH 

Fig. 10. Variation of oonthly aean rainfall (■■) averaged between 30*11 and 

30* S (fron Jaeger, 1976), and Intenalty of aubtroplcal hlgha a«aaured 
by nuaber of (4* lat. x 5* long.) grid polnta for which aea level 
presaure la greater than 1020 nb (froa Godbole and Shukla, 1981). 
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Fig. 11. Stationary variance of January scan geopotential height at 500 ab 

■ f^BulAted by two versions of the GLAS cllaute aodel, and observations 


